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I.   INTRODUCTION 


A.   EACKC-FCUJiE 


As  cf  1S74,  the  United  States  contained  only  6%  cf  the 
world  population  but  was  using  33%  of  the  energy  consumed 
each  year  according  to  Ref.  1.  Approximately  one  third  of 
the  oil  used  in  the  U.S.  is  imported.  This  heavy  dependence 
en  foreign  scurces  fcr  energy  coupled  with  the  exhaustion  of 
dcuestic  fessile  fuels,  particularly  oil,  in  the  foreseeable 
future  is  causing  the  U.S.  military  and  civilian  sectcrs  of 
the  economy  tc  search  for  energy  sources  that  wculd  be 
inexhaustible  and  independent  of  foreign  ccntrol. 
Increasing  fuel  costs  and  environmental  prchlems  are  ether 
forces  pushing  the  search  for  alternatives.  Nuclear  power, 
once  considered  the  cure-all,  is  facing  severe  problems. 
Ccsts  cf  censtruction,  operation,  and  fuel  have  risen 
dramatically.  Questions  concerning  the  safety  of  operation 
of  the  nuclear  plants  and  of  the  storage  cf  nuclear  vastes 
ha?e  beccne  pclitical  issues  that  have  slowed  the 
exploitation  cf  nuclear  power.  The  present  nuclear  reactors 
use  a  fuel  that,  like  petroleum,  has  a  limited  availakility 
and  whose  price  is  increasing  rapidly.  The  breeder  reactor, 
which  would  create  more  fuel  than  it  burns,  and  the  fusion 
reactor,  which  would  use  the  hydrogen  foucd  in  water,  are 
still  a  leng  way  frcm  producing  power  for  a  world  whose 
deoands  fcr  energy  increase  each  year. 

Many  people  advocate  the  exploitation  cf  the  sc-called 
"free"   enercry  scurces.   These  scurces  are:   (1)  the  kinetic 
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energy  cf  moving  fluids,  such  as  winds  and  ccean  currents; 
(2)  the  pctential  energy  of  tides  and  rivers;  (3)  the  heat 
generated  inside  the  earth;  (4)  the  direct  conversion  of 
solar  energj  into  electricity  and  heat;  and  (5)  the  use  of 
sun  warmed  water  in  conjunction  with  a  source  of  colder 
water  tc  provide  the  temperature  differential  to  run  a  vapor 
power  cycle.  None  cf  these  sources  are  truly  "free".  fihat 
is  happening  is  the  trading  off  cf  the  transportation, 
processing,  and  environmental  costs  of  a  concentrated  energy 
source  for  the  capital,  operating,  and  social  costs  of 
converting  a  very  diffuse  source  into  a  mere  concentrated 
form  such  as  electricity. 

Ihe  earth  can  te  thought  of  as  a  giant  heat  engine, 
absorbing  erergy  frcm  the  sun  and  reflecting  energy  as 
thermal  radiaticn  rack  to  space  as  pictured  in  Bef.  2. 
Since  the  pcles  receive  less  energy  than  the  equator,  the 
atmosphere  and  ocean  attempt  tc  distribute  the  energy  more 
evenly  ever  the  earth.  The  air  is  heated  ty  absorption  of 
sunlight  and  ty  contact  with  the  surface  of  the  earth.  The 
heated  air  rises  thus  setting  up  surface  wind  currents  as 
cooler  air  tries  tc  replace  the  rising  air.  The  shear 
stress  created  by  the  relative  motion  between  the  air  and 
the  surface  water  causes  the  water  to  move  thus  creating  the 
surface  water  currents.  Many  currents  in  the  oceans, 
regardless  cf  depth,  originate  with  the  wind  shear  stress  at 
the  air-surface  interface.  The  circulaticn  patterns  cf  the 
oceans  beccme  extremely  complex  due  to  the  influence  cf  the 
rotation  cf  the  earth  and  the  shapes  of  the  ccean  basins. 

In  general  the  surface  waters  near  the  equator  are 
warned  by  the  sun  and  flow  towards  the  poles  giving  up 
energy  by  evaporation,  and  by  convection  tc  the  air  ard  by 
radiaticn  rack  to  space.  As  these  currents  cool,  seme  of 
the  water  becomes  dense  enough  tc  sink  and  then  flows  back 
to  the  eguatcr  alcng  the  bottom. 
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The  relationship  between  salinity  and  temperature  and 
the  density  cf  seawater  is  such  as  tc  create  a  stable 
stratif icaticn  of  the  ocean  in  many  areas.  A  salinity 
decrease  cr  a  temperature  increase  reduces  the  density.  The 
total  effect  is  to  create  a  warm  surface  current  overlaying 
the  cold  water  returning  from  the  poles.  In  the  tropical 
waters  the  surface  temperature  may  exceed  80°F  (26.7°C) 
while  3,00C  to  6,000  ft.  (914  to  1828  m)  dewr  the 
temperature  nay  be  about  40°f  <4.4°C).  Ihis  temperature 
difference  can  te  used  to  rur  a  roan-made  heat  engine,  like 
the  Ocean  Thermal  Energy  Converion  (OTEC)  power  plant 
ccrcept. 

The  proper  choice  of  a  site  is  critical  to  the  cost  of 
such  a  heat  engine.  Since  the  temperature  difference  is  low 
ccupared  to  conventional  thermal  power  generating  methods 
which  typically  use  temperature  differentials  exceeding 
50C°F  (26C°C),  the  effect  of  a  1°F  (0.56°C)  loss  in  the  AT 
is  much  nore  costly  at  the  low  thermal  differentials,  which 
are  usually  less  than  50°F  (27.8°C)  in  the  ccean. 
Therefore,  the  vertical  temperature  profile  of  the  cceans 
will  ccntrcl  the  location  of  an  OTEC  power  plant.  The 
surface  temperatures  vary  with  location  and  with  time.  At 
mest  places  in  the  ocean  there  is  a  seasonal  change  of  the 
surface  later  temperature  of  5  to  8°F  (2.8  to  4.4°C)  except 
in  the  trcpical  areas  that  remain  relatively  constant 
according  to  Bef.  3.  However,  all  regions  of  the  oceans  are 
affected  ty  wind-yenerated  mixing  of  the  surface  water  that 
changes  the  surface  temperature  by  a  degree  or  so  and  a 
daily  cycle  cf  a  few  tenths  cf  a  degree.  Some  cf  the 
possible  0 . £ .  near-shore  sites  being  considered  are  the 
waters  arcund  the  Kawaian  Islands  and  the  Gulf  Streai  off 
the  southeastern  U.S. 

In  the  past,  the  engineering  analysis  and  the  cost 
estimation   have   teen   conducted  as   if  they  were  separate 
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functions.  In  a  system  with  a  snail  temperature 
differential,  such  as  this  one,  the  engineering  design  and 
economic  analysis  should  be  linked  together  as  a  single 
procedure  because  of  the  capital  intensive  nature  cf  the 
ocean  thermal  power  plant.  For  example,  the  heat  transfer 
egiaticn  (Q  =  UA  AT)  is  coupled  to  the  cost  equation  for 
heat  exchangers  ($  =  K  An )  through  the  surface  area  of  the 
tutes.  Eecause  of  the  small  AT  involved  the  area  must  be 
very  large  tc  supply  the  same  Q.  The  other  components  are 
all  linked  together  by  the  engineering  and  cost  equations  so 
that  the  picblem  cannot  logically  be  troken  intc  an 
eguipment  by  equipment  optimization.  The  present  analysis 
cocsiders  a  closed  vapor  power  system  operating  cu  the 
thermal  pctertial  available  from  the  ocean.  The  economic 
optimization  cf  the  system  is  carried  out  to  shew  what 
information  can  be  gotten,  what  conclusions  can  be  reached 
from  such  analysis  and  where  research  effort  shculd  be 
expended  tc  improve  the  design  of  the  system.  This  study 
dees  not  attempt  tc  predict  the  cost  of  such  a  system  or  tc 
make  specific  recommendations  about  the  parameters. 

A  schematic  diagram  of  the  system  is  shewn  in  figure  1. 
The  warn  surface  water,  the  Gulf  Stream  for  example,  is 
pumped  through  a  heat  exchanger.  There  seme  energy  is 
transferred  by  boiling  the  working  fluid  such  as  ammonia. 
The  working  fluid  is  piped  to  a  turbine  where  some  of  the 
energy  is  ccryerted  to  mechanical  energy  then  to  electcical 
energy  in  a  generator.  The  unavailable  portion  of  the 
absorbed  energy  is  rejected  to  the  environment  in  the 
condenser  and  the  working  fluid  is  returned  to  the  liguid 
state.  The  lew  temperature  in  the  condenser  is  maintained  by 
pumping  ccld  ocean  water  up  from  the  depths.  Finally  the 
working  fluid  is  pumped  back  to  the  boiler. 

The  Ccear  Thermal  Energy  Conversion  program  has  received 
much    attention    in   the   last   several   years   from   many 
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researchers  funded  at  first  primarily  by  the  National 
Science  Foundation  (NSF) .  Since  the  establishment  cf  the 
Energy  Feseaich  and  Development  Administration  (ERDA)  ,  new 
research  contracts  are  funded  mainly  by  this  agency.  The 
l.S.  Navy  is  alsc  funding  some  research  in  the  OTEC  area. 
Ref.  4,  5  and  6  give  a  geed  summary  of  the  work  that  has 
been  done.  fcbile  much  of  the  work  has  concentrated  or  the 
heat  exekancers,  the  other  components  and  the  systen  as  a 
whele  are  receiving  attention.  Some  researchers  are 
concentrating  on  the  environmental  and  legal  guestiens  of 
such  systems,  ether  teams  are  considering  hew  to  use  the 
energy  supplied  from  an  OTEC  power  plant. 

In  any  design,  trade-offs  must  be  made.  Capital  cost 
can  be  traded  for  operational  costs.  Pumping  power  can  be 
traded  fcr  heat  exchanger  size.  Trade-cffs  can  involve 
items  that  have  a  common  measure  such  as  dollars;  seme  ether 
trade-offs  are  very  hard  to  measure  in  common  terms  such  as 
social  and  environmental  costs.  The  trade-cffs  may  involve 
crly  one  cempenent  such  as  a  pump  or  a  turbine,  but  most 
often  the  tiade-effs  affect  more  than  one  component  or 
possibly  the  entire  system. 

The  trade-eff  process  should  continue  until  an  cctimum 
design  is  obtained.  From  Ref.  7,  "The  aim  cf  optimization 
is  the  selection,  cut  of  the  multiplicity  of  potential 
solutions,  of  that  solution  which  is  the  best  with  respect 
tc  some  well-defined  criterion."  In  order  for  the  plant 
design  tc  D€  optimum,  the  components  must  be  optimized 
within  the  context  cf  the  entire  systen.  Some  cf  the 
parameters  fcr  a  piece  of  eguipment  may  have  no  direct 
effect  en  the  optimum  design  of  another  portion  cf  the 
system,  but  seme  parameters  will  have  significant  effects. 
For  examrle,  a  heat  exchanger  is  made  up  cf  a  large  number 
of  tubes  cf  a  given  diameter  and  length  and  is  designed  to 
have   a   certain  heat  transfer  rate  for  a  given  flow  rate  of 
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water  pumped  through  the  tubes.  The  diameters,  length  and 
the  numter  cf  tubes  may  be  optimized  for  lowest  cost  cf  the 
heat  exchanger.  This  approach  leaves  out  the  effect  that 
the  tube  diameter  and  the  tube  length  have  en  the  pressure 
drcp  across  the  heat  exchanger.  Minimizing  the  cost  cf  the 
heat  exckaEger  irvclves  minimizing  the  surface  area.  This 
may  lead  tc  increasing  the  pressure  drcp  acrcss  the 
exchanger.  Ihe  higher  the  pressure  drop  the  more  powerful 
and  the  irere  expensive  the  pump  must  be  to  maintain  the  flow 
velccity.  Ibis  is  just  one  example  cf  the  hundreds  of 
trade-offs  pcssible. 

Onf orturately,  what  happens  in  many  cases  is  that  each 
grcup  cf  experts  working  on  a  particular  component  optimize 
that  iten  arcund  some  given  conditions  that  may  have  teen 
set  ty  teams  working  en  some  other  section  cf  the  system. 
Little  analysis  is  made  of  how  the  entire  system  responds  to 
the  coupling  constraints.  All  that  the  puirp  manufacturer 
wants  is  that  the  buyer  tell  him  what  the  flow  rate,  head 
and  service  conditions  are  and  he  will  guote  a  price. 
Before  the  buyer  can  decide  on  the  flow  rate  and  head  he 
should  first  know  how  the  ccst  cf  the  pump  varies  as  flow 
rate  and  head  vary  so  that  those  two  parameters  can  be 
optimized  in  the  ccntext  of  the  system. 

The  prctlem  then  becomes  one  of  acguiring  sufficient 
information  en  hew  costs  vary  as  certain  parameters  change. 
Each  manufacturer  is  able  to  make  estimates  cf  costs  fcr  his 
picduct  tut  re  is  reluctant  to  devulge  his  informaticr  to 
buyers  cr  cempetitors.  This  makes  constructing  cost  cur /es 
a  hazardcus  prccess  at  best,  since  the  data  must  be  gathe— ?d 
by  ether  nethcds  that  may  cause  large  unexplained  variances. 

Sometimes  the  optimization  analysis  is  made  with 
assumpticES  that  make  the  answers  cf  guestionable  value. 
Fcr  example,  sometimes  the  heat  transfer  coefficient,  U,   is 
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assumed  fixed  in  the  equation  q  =  UA  AT.  But,  as  shewn  in 
Ref.  8,  for  toilinq  the  heat  transfer  coefficient  is  very 
sensitive  tc  the  temperature  difference  with  the  result  that 
q  =  UA'Al3-33.  The  The  accuracy  of  the  engineering  eguations 
should  te  kept  in  mind.  Most  of  the  eguations  are 
correlations  tc  fit  experimantal  data.  The  correlation  was 
developed  for  seme  particular  set  of  data  and  may  disagree 
signif icartly  when  compared  to  data  taken  by  ether 
experimenters.  The  limitations  of  the  various  relations 
must  be  kept  in  mind  when  reviewing  the  results.  for  a 
system  the  size  of  the  OTEC  plant,  a  pilot  plant  should  be 
constructed  in  order  to  prove  out  the  answers  given  bj  the 
calculations  before  expensive  mistakes  are  made. 

As  a  design  progresses,  the  analysis  begins  with  a 
general  picture  containing  many  simplifying  assumptions  and 
proceeds  foreward  with  more  details  added  at  each  stage. 
This  report  is  the  second  one  written  at  the  Naval 
Pcstraduate  fchcel  on  the  subject  of  OTEC.  Commander  Jurman 
Sheppard,  C£K,  in  his  analysis,  Ref.  9,  attempted  tc  show 
what  kinds  cf  information  could  be  gained  from  a  ccnbined 
thermal  eccccnic  model.  He  advocated  the  use  cf  zone 
analysis,  as  developed  by  Wismer  in  Ref.  10,  fcr  large 
complex  systens.  In  Wismer* s  approach  to  zone  analysis,  the 
system  is  brcken  dewn  into  zones  containing  one  or  more 
cempenents.  She  zones  are  connected  by  linking  variables. 
Two  metheds  cf  optimization  are  possible.  For  the  first 
method,  in  each  zone,  the  zone  parameters  are  optimized 
while  the  linking  variables  are  held  fixed.  Then,  the  zone 
parameters  are  held  fixed  while  the  linking  variables  are 
optimized.  lie  process  iterates  between  optimizing  the  zone 
parameters  and  the  linking  variables  until  convergence 
criteria  are  satisfied.  In  the  second  method  the  zenes  are 
cut  off  ticn  each  other  and  each  linking  variable  is 
separated  into  two  variables,  one  on  each  side  cf  the 
cutting   flare.   Tie  optimization  proceeds  until  the  linking 
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variables  cr.  each  side  of  the  cuts  are  egual. 

Sheppard's  system  did  not  produce  vapor  in  a  boiler  tut, 
instead,  maintained  the  working  fluid  in  the  liguid  state  in 
the  wain  seawatei  heat  exchanger  and  vaporized  the  working 
fluid,  aancnia  in  his  model,  by  expanding  it  through  a 
"black  fccx"  called  a  vaporizer.  His  optimization  analysis 
was  confined  to  the  zone  containing  the  warm  seawater 
exchanger,  the  feed  pump,  and  the  warm  water  circulating 
pump. 


E.   CEJEC!IIVIS 


The  present  analysis  improves  the  previous  model 
analysis  ky  replacing  the  liguid-to-liguid  heat  exchanger 
and  the  vaporizer  with  a  boiler.  The  boiler  model  is  to  be 
of  sufficier.t  detail  so  that  a  realistic  desigr.  is 
represented.  The  ether  major  components  of  the  system,  such 
as  the  turbire,  the  condenser,  and  the  condenser  circulating 
pumps,  are  included.  A  nonlinear  programming  technique  is 
tsec  to  perform  the  thermal-economic  optimization  of  the 
completed  system.  The  research  reported  here  is  intended  to 
show  the  usefulness  of  optimization  analysis  in  the  design 
of  a  lew  thermal  difference  system  and  to  demonstrate  some 
of  the  problems  that  must  be  considered  in  this  type  of 
analysis . 
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II.   THE  MODEL 


A.   GESE5II  lESCBItTION 


Vapcr  pcher  generating  cycles  usually  consist  of  four 
processes.  Energy  is  absorbed  from  a  high  temp eiature 
source  and  is  used  to  vaporize  a  working  fluid  in  a  toiler. 
Next,  work  is  extracted  by  expanding  the  vapcr  in  an  engine, 
exbaustinc  tie  vapcr  at  a  lower  temperature  and  pressure. 
The  working  fluid  is  then  condensed  back  to  the  liguid  state 
by  rejecting  energy  to  the  lcwer  temperature  energy  sink  in 
the  third  process.  Finally,  the  liquid  is  returned  to  the 
higher  pressure  of  the  boiler  by  a  pump.  In  the  following 
analysis,  the  energy  source  is  a  current  of  ocean  water 
heated  by  tbe  sun  and  the  energy  sink  is  assumed  to  be 
cclder  ocean  water  from  deep  currents.  The  hcrking  fluid  is 
assumed  tc  he  propane  although  many  other  fluids  are 
possible,  such  as  ammonia  or  the  Freons. 

Six  ccirpcnents  are  modeled  in  the  following  analysis. 
The  hot  side  heat  exchangers,  the  boilers,  vaporize  the 
working  fluid  which  is  expanded  through  a  turbine  exhausting 
to  a  ccndenser  that  rejects  the  unused  heat  to  the  cold 
sirk.  Three  pumps  are  included,  a  feed  pu&p  returns  the 
wcrking  fluid  to  the  boiler  and  two  pumps  circulate  the  warm 
and  cold  waters  through  the  boiler  and  the  ccndenser, 
respectively.  Ihe  arrangement  of  the  system  is  shown  in  Fig. 
2.  The  ncdes  are  numbered  to  correspond  with  the  fluid 
state  points  shewn  later,  in  Fig.  3. 
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The  analytical  model  consists  of  the  eguations 
describing  the  system  and  a  set  cf  parameters  to  be 
optimized  sc  as  tc  minimize  the  capital  cost  of  the  plant 
subject  tc  certain  constraints.  The  objective  function,  the 
functicn  tc  be  minimized,  is  the  sum  of  the  capital  costs  of 
the  pumps,  the  boiler  and  the  condenser.  Although  the 
turbine  is  a  significant  capital  cost  and  would  affect  the 
design  cf  the  system,  it  is  not  included  in  this  nodel. 
There  are  twc  classes  cf  constraints.  The  first  class  is  the 
explicit  ccnstraints  that  specify  the  upper  and  lower  bounds 
on  the  parameters.  The  other  class  of  constraints  is  the 
implicit  ccnstraints.  These  are  the  engineering  eguations 
that  describe  the  operation  of  the  power  plant.  For 
example,  a  ccnstraint  might  say  that  the  teaperature  inside 
the  boiler  rorst  be  less  than  the  temperature  of  the  warm 
seawater  because  the  eguations  fail  if  the  reguirement  is 
net  satisfied.  Another  example  would  be  a  power  balance 
that  reguires  that  the  sum  of  the  energy  flews  into  and  out 
of  the  systen  be  zero.  Additionally,  they  proscribe  the 
regicn  that  the  cptimizaticn  rcutine  may  search  within  for 
the  miniffum  cf  the  problem.  The  implicit  ccnstraints  are  of 
twc  types,  eguality  and  ineguality  constraints.  Ineguality 
ccnstraints  set  one-sided  bounds  on  the  feasible  region. 
Seme  of  the  inequality  constraints  are  restrictions  placed 
on  the  maximum  or  ninimum  values  of  the  parameters  due  to 
judgements  cade  en  the  basis  of  practical  engineering 
consider atic rs  .  The  equality  constraints  are  relationships 
that  must  be  strictly  satisfied  at  the  optimimum.  The  system 
of  eguatiens  describing  the  model  fit  intc  the  following 
f  cut : 
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f  (X) 


subject  to:   g  (X)  >  0,  i=  1,2, .  . . .  ,m 

i 

h  (X)  =  0,  i  =  m+ 1 , .  . . . ,m+n 
i 


where  X  is  the  vector  of  the  parameters  tc  be  optiiiized, 
f (X)  is  the  capital  cost  function,  g^  (X)  are  the  inequality 
constraints  and  hj  (X)  are  the  equality  constraints.  The 
parameters  are  the  dimensions  of  the  toilers  and  the 
condensers,  the  boiling  and  condensing  pressures,  the  sea 
water  velocities,  and  the  mass  flew  rate  of  the  working 
fluid  through  the  system. 

The  tasis  for  the  economic  framework  must  be  set  cut  at 
the  start.  £11  measurable  costs  are  relevant  to  the  final 
design  of  the  power  plant,  but,  for  some  types  of  analysis 
and  design,  irany  ccsts  can  be  ignored.  The  life  cycle  cost 
of  the  CTEC  must  be  considered  before  the  decision  is  Bade 
tc  commit  the  U.S.  and  the  rest  of  the  world  tc  the 
censeguences  cf  choosing  this  energy  generation  scheme  as 
even  a  partial  solution  to  the  energy  problem.  The  life 
cycle  ccst  consists  of  the  capital  costs,  the  financing 
costs,  and  the  operating  costs.  The  trade-off  of  operating 
ccst  and  capital  cost  is  not  considered  in  this  research, 
because  the  operating  cost  is,  for  the  most  part, 
independent  cf  the  overalL  dimensional  characteristics  of 
the  compenents.  Operating  costs  are  dominated  bj  the 
structure,  site  location,  and  maintenance  requirenents. 
Such  things  as  the  kind  of  instrumentation,  controls, 
bearings,  and  structural  materials  have  mere  effect  en  the 
operational  ccsts  than  the  length  of  the  bciler  or  the 
diameter  cf  the  toiler  shell.  Therefore,  enly  the  capital 
ccsts  of  the  major  components,  except  the  turbine,  are 
considered . 
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To  fc€  iseful  for  model  analysis  cost  data  must  be 
transformed  into  a  cost  estimating  relationship  (CER)  that 
shews  he*  cost  varies  with  some  variable  or  some  comhication 
of  variables.  The  costs  must  be  demonstrated  to  be  highly 
correlated  with  the  variable  chosen.  Since  the  cost 
relation  is  enly  a  correlation,  the  bounds  on  tne  accuracy 
of  the  equation  shculd  be  specified  by  those  who  develop  it. 
Since  the  relation  is  arrived  at  from  a  finite  body  cf  data, 
the  CER  is  directly  applicable  only  in  the  range  of  the  data 
collected.  Extrapolation  cf  the  CER  outside  the  rar.ge  is 
cargercus  because  there  is  nc  data  in  that  region  tc  support 
the  assumed  curve.  However,  the  value  of  CER's  is  the  use 
of  them  tc  enable  the  estimator  to  make  a  better  analysis  of 
a  new  design  than  he  could  if  the  information  is  not 
available.  The  relevant  range  of  the  CER  must  be  kept  in 
nind  wher  analyzing  the  solutions  from  the  standpoint  cf  the 
reliability  cf  the  answers.  Another  factcr  tc  consider  when 
using  cost  estimating  relationships  is  whether  the  equipment 
the  data  was  taken  from  is  comparable  to  the  systen  the 
analyst  is  considering.  If  the  design,  the  technology,  or 
the  application  of  the  proposed  equipment  is  significantly 
different  ficm  that  from  which  the  data  is  taken,  the  CER 
may  give  unreliable  answers. 

This  aralysis  begins  with  the  development  of  the 
engineering  eguatiens  that  constrain  the  system.  This  is 
followed  by  the  development  of  a  cost  framework  based  on 
capital  costs  as  the  function  to  be  minimized.  The 
engineering  and  economic  relations  are  linked  together  in 
the  objective  function  and  the  constraints.  A  description 
of  the  algorithm  that  performs  the  optimization  is  included. 


B.   CYCLE  AtiMSSIS 
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The  tuilcing  of  the  model  begins  with  the  analysis  of 
th€  ther ncdynamic  cycle.  Figure  3  is  a  temperature-entropy 
diagram  of  the  working  fluid  on  which  the  cycle  is  shown. 
These  numbered  ncdes  represent  the  thermodynamic  state 
pcints  at  that  location  in  the  system  and  correspond  tc  the 
numbers  ir  Big.  2.  The  boiling  and  condensing  processes  are 
assumed  tc  be  reversible  processes.  In  the  pump  arc  the 
turbine  the  reversible  process  is  indicated  by  the  subscript 
"s"  on  the  state  pcint  number.  In  the  actual  system  there 
are  irreversible  losses  inside  each  component  that  are  left 
out  of  the  analytical  model  or  lumped  in  with  ether  losses. 
The  pressure  drops  due  to  frictional  losses  in  the  toiler 
and  the  cendenser  are  neglected. 

State  pcint  1  assumes  that  the  liquid  leaving  the  toiler 
is  saturated.  State  points  2s  and  2  represent  the 
conditions  leaving  the  feed  pump.  If  there  were  no  lesses  in 
the  pump,  the  fluid  would  be  at  2s.  The  fluid  enters  the 
bciler  tc  he  heated,  first,  to  saturation  conditions  at 
state  point  3  then  boiled  at  essentially  constant  pressure. 
State  pcint  4  shows  the  vapor  leaving  the  bciler  as  a 
saturated  vapor.  Ihe  lines  from  point  4  to  5s  and  5 
represent  the  expansion  through  the  turbine.  The  path  tc 
state  pcint  5  accounts  for  the  losses  due  tc  frictior  and 
ether  irreversitle  processes  in  the  turbine.  The  remaining 
heat  is  then  rejected  from  point  5  to  point  1,  again  assumed 
to  he  a  constant  pressure  process. 

Ihe  theoretical  heat  transfer  rate  in  the  boiler  and  the 
cendenser  is  found  by  subtracting  the  enthalpy  at  state 
pcints  2  and  1  from  the  enthalpy  of  state  pcints  4  ar.d  5, 
respectively,  and  multiplying  by  the  total  mass  flow  rate  of 
the  working  fluid. 

Qab  =  M^(H4  -  Hz)  (1) 
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FIGURE  3.  THERMODYNAMIC  CYCLE 
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Qrcj  =  MW|(H5  "  H,) 
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The  enthalpy  at  state  point  5  is  determined  by  first  finding 
the  moisture  ccntent  of  the  fluid  exiting  the  turbine. 
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The  enthalfy  at  point  5s  is  determined  from 


<3) 


H5s  ~  H1  +  i  H*cj 


I**) 


This  allows  E5  to  te  found  from 


'tur  ~ 


H4  -  H5 
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where   the   efficiency   of  the  turbine  is  assumed  kncwr.  The 
power  out  of  the  turbine  is  now  known. 


Ptur  -  Mu/-f  <HH  "  H5) 


(6) 


For  the  feed  pump  the  change  in  enthalpy  is  accounted  for  by 
the  change  ir  pressure  assuming  negligible  rise  ir  the 
temperature.  The  enthalpy  of  point  2  is  determined  free 
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H2  -  H1 


(7) 


where  the  efficiency  is  assumed  known.  The  power  reguired  to 
drive  the  {unp  is  determined  from 
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Up     =  Mw-MH2  "  H1) 


<3) 


Tc  gain  a  better  feeling  fcr  the  low  thermal  potential 
systems  end  tc  irak€  comparisons  with  other  systems  several 
efficiencies  are  worthwhile  calculating.  The  first  is  the 
Carnct  efficiency  which  defines  the  maximum  efficiency 
possible  fci  a  system  operating  between  two  temperature 
Units  ir  decrees  absolute,  °H  (°K). 
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The  next  efficiency  calculation  takes  into  account  the 
theoretical  tnermcdynamic  cycle  and  is  the  ratio-of  the  net 
pofcer  cut  cf  the  cycle  to  the  heat  taken  up  ty  the  cycle 
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This  acccints  for  the  thermodynamic  inefficiencies  inherent 
in  the  cycle.  Finally,  the  system  efficiency  takes  into 
account  the  ether  internal  consumers  cf  power:  tho 
circulaticn  pumps,  the  pipe  losses,  and  any  other  lesses, 
such  as  the  antifouling  equipment.  This  calculation 
acccunts  fcr  all  of  the  mechanical  and  electrical 
ine i f iciercies  of  the  system. 
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For  example,  if  the  hot  water  is  at   75°F   (23.9°C)   and 
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the  ccld  water  is  assumed  at  40°F  (4.4°C)  the  Carnot 
efficiency  is  6.5%.  Per  one  typical  set  of  assumptions,  as 
shewn  in  Table  5,  where  the  turbine  and  the  pumrs  are 
assumed  tc  be  65%  efficient,  the  cycle  efficiency  is  2.6? 
and   the  system  efficiency  is  2.4%   for  one  feasible  design. 


C.   ENGIKEEB1NG  FRAMEWOBK 


This  secticn  describes  the  design  of  the  components  that 
make  up  the  system.  The  level  of  detail  of  the  analysis 
prcvides  a  physically  realistic  model  and  clcses  the  lcop. 
Ike  pipe  friction  losses  on  the  cold  and  hct  water  inlets  as 
well  as  tfce  fipe  losses  in  the  working  fluid  side  are 
neglected.  ihe  turbine  is  included  only  as  an  energy 
transducer  and  no  attempt  is  made  to  obtain  costs,  physical 
diirensicES,  cr  fluid  flow  characteristics.  The  efficiencies 
of  the  turbine  and  the  pumps  are  assumed  as  known  constants. 
The  heat  arscrbed  in  heating  the  working  fluid  tc  the 
saturated  state  is  a  small  portion  of  the  tctal  energy  taken 
uc  in  the  bciler;  therefore,  the  sensible  heating  cf  the 
fluid  is  neglected. 

1 .   lie  Eciler 

The  tciler  consists  cf  a  number  cf  tube-in-shell 
heat  exchangers  connected  in  parallel  as  shewn  in  Fie.  2. 
Figure  4  shows  the  internal  arrangement  of  cne  shell.  The 
tube  bundle  is  horizontal  and  boiling  occurs  on  the  outside 
surface  cf  the  tubes.  The  tube  bundle  is  submerged  in  the 
working  fluic.  The  working  fluid  is  maintained  at  a  le^el  to 
just  cover  the  top  rcw  of  tubes.  The  warm  seawater  is  tumped 
through  the  inside  of  the  tubes. 
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FIGURE  4.  BOILER  CROSS  SECTION 
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lhe  clearance  between  the  tubes  is  fixed  at  0.5 
inches  (1.27  cm)  to  allow  for  flaring  of  the  tube  ends  or 
ether  space  and  strength  reguiiements.  Similarly,  the  tube 
thickness  is  fixed  at  0.028  inches  (0.071  cm),  since  no 
strength  requirements  due  to  differences  between  internal 
and  external  pressures  on  the  tubes  or  material  loss  due  to 
ccrrosion  cf  the  tube  surfaces  are  accounted  for. 

At  the  tcp  of  the  tube  bank  the  leaving  vapor 
certains  a  large  amcunt  of  liquid  droplets  which  must  be 
eliminated  eitfcer  in  the  boiler  or  by  an  external  separator. 
In  the  mcdel  the  separation  takes  place  in  the  bciler  by 
allowing  sufficient  space  above  the  tube  bank.  Ealen  and 
Small,  Ref.  11,  presented  a  method  for  finding  the  necessary 
separaticc  volume  that  considers  the  effects  of  surface 
tension,  density  of  the  vapor  and  the  liquid,  and  the  mass 
flew  cf  tbe  wcrking  fluid.  First  the  allowable  vapor  lead  is 
found  frcn 


VL  = 


2290^  [- 


0.5 


86x10-5  (^  -  ps) 


(12) 


vapcr  load  is  an  expression  for  the  mass  flow  rate  cf  vapor 
through  a  unit  volume.  The  expression  for  the  maximum  vapor 
lead  limits  the  flew  rate  so  that  the  moisture  droplets  are 
able  to  settle  out  of  the  vapor  before  it  leaves  the  bciler. 
Then  the  reguired  segmental  area  is  found  frcm 


SA  = 


«W 


VL  L 


(13) 


The  segmental  area  is  the  cross  sectional  area  from  the  top 
of  the  tute  bank  to  the  top  of  the  shell.  Kncwing  the 
segmental  area  the  height  of  the  tube  bank   above   the   main 
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diameter  is  fcund  from 


SA  =  C.25  E£   ARCOS 


m- 


2  2  0.  5 

hs(0.25  Is      -  fas  )        (14) 


Falen  and  Snail  also  present  a  much  simpler  method  of 
determining  the  total  height  of  the  tube  bank.  The  height 
of  the  separation  volume  is  set  at  U0%  cf  shell  diameter 
which  results  in  nhs"  being  equal  to  0.1  Ds.  Either  nethod 
car  be  used  alone  or  the  first  method  can  be  usee  as  a 
censtraint  on  the  second.  For  this  protlem  the  segnental 
area  is  fcund  by  bcth  methods  with  the  first  used  as  a  check 
en  the  values  given  by  the  second. 

Since  the  heat  flux  calculation  is  based  on  a  single 
tube,  the  nunber  of  tubes  per  shell  must  be  determined  in 
order  tc  find  the  total  energy  absorbed  by  the  heat 
exchanger.  7he  tubes  are  assumed  to  be  laid  out  in  a  60 
degree  triangular  pattern  for  maximum  compactness.  The 
tures  are  lecated  at  the  vertices  of  an  equilateral 
triangle.  iwc  cross  sectional  areas  are  calculated.  The 
area  of  a  hexagon  whose  minor  radius  is  egual  to  cne  half 
the  pitch  is  determined  from 


=  6 


©■ 


TAN  30° 


(15) 


Then  the  effective  cross  sectional  area  of  the  tube  bank  is 
fcund  by  subtracting  the  clearance  area,  SA,  from  the  cross 
sectional  area  of  the  shell  after  correcting  the  shell 
diameter  for  a  clearance  space  to  allow  for  free  flow 
channels  both  within  the  tube  bundle  and  arcund  the  ottside 
of  the  tuhe  bundle.  The  effective  cross  sectional  area  of 
the   tube   bark  is  divided  by  the  area  of  the  hexagen,  found 
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afccve,  tc  determine  the  number  of  tubes  in  a  shell. 

Ire  flew  of  the  liquid  and  vapor  through  the  tube 
fcar.k  is  assuned  to  be  turbulent.  The  incoming  sufccccled 
working  fluid  from  the  feed  pump  uixes  with  the 
recirculating  fluid  returning  frcm  higher  up  in  the  tube 
bank.  The  circulation  of  the  fluid  is  induced  ty  the 
boiling  process.  This  mixture  enters  the  tube  tank  and 
begins  the  bciling  process.  Near  the  bottom  of  the  tube 
bark  the  heat  transfer  goes  mostly  to  bring  the  fluid  up  to 
the  boiling  temperature.  After  the  fluid  reaches  saturated 
liguid  cenditiens,  stable  boiling  begins.  The  vapor  bubbles 
rising  thrcugh  the  surrounding  liguid  increase  the 
turbulence  aicund  the  tubes.  As  the  liquid  continues  up 
thrcugh  the  tube  bank,  the  local  boiling  temperature  is 
decreasing  due  to  the  decreasing  hydrostatic  head.  This 
decrease  in  the  boiling  temperature  increases  the 
ten perat ure  difference  and  increases  the  heat  fie*  and 
increases  the  bubble  generation  rate.  Tne  vapor  generation 
may  become  sc  rapid  that  a  condition  called  vapor  blanketing 
may  occur  in  the  interior  of  the  tube  bank  where  the  liquid 
cannot  flow  inward  fast  enough  to  displace  the  vapor  flowing 
through  the  region.  This  is  similar  to  burnout  in 
cenvective  bciling  in  that  the  tubes  dry  out  and  the  heat 
transfer  rate  decreases.  This  effect  becomes  proncunced 
enly  near  the  critical  heat  flux  according  to  Starczewski, 
in  Ref.  6. 

All  cf  the  heat  and  nass  flow  variations  in  the  tube 
bundle  are  averaged  by  performing  the  calculations  for  a 
single  tube  located  at  the  mid  height  of  the  tube  bank. 
This  ture  is  assumed  to  represent  the  average  of  all  the 
tubes  in  the  tube  bank.  The  boiling  situation  is  assuned  to 
be  saturated  pool  boiling.  Quoting  frcm  Collier,  Ref  12, 
"pool  boiling  is  defined  as  boiling  from  a  heated  surface 
submerged   in  a  large  volume  of  stagnant  liquid".   Since  the 
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fluid  vcluire  seen  ty  the  tube  is  neither  large  nor  stacnant, 
a  ccrrecticc  is  reeded  for  the  effect  cf  being  in  a  tube 
bundle.  Saturated  boiling  is  defined  as  boiling  at  a 
constant  temperature  and  pressure  as  the  liguid  changes  to  a 
vapor.  Ihis  is  probably  not  strictly  true  because  the 
boiling  temperature  is  controlled  by  the  local  pressure  at 
each  tute,  and  the  fluid  is  not  reaching  a  stable 
temperature  as  it  rises  in  the  bundle.  The  situation  is 
probably  closer  to  convective  boiling  since  the  fluid  is 
boiling  in  a  confined  channel  between  the  tubes.  An 
additional  effect  in  the  boiler  is  the  fact  that  the 
hydrostatic  head  at  the  bottom  of  the  tute  bundle  causes  the 
boiling  tc  teke  place  at  a  higher  pressure  at  the  tottom 
than  the  toiling  at  the  top  of  the  bundle  where  the  pressure 
is  egual  tc  the  boiler  exit  pressure.  This  effect  is 
accounted  for  by  using  the  average  head  of  the  bundle  in  the 
heat  trarsfei  egcatiens. 

Ihe  heat  transfer  analysis  of  the  boiler  begins  with 
the  tasic  eauation 


Q  =  UA  AT^ 


(16) 


where  U  is  the  overall  heat  transfer  coefficient,  A  is  the 
area  thrcugfc  which  the  heat  is  transferred  and  ATm  is  a 
properly  defined  mean  temperature  difference.  In  this 
paper,  the  product  UA  is  treated  as  a  single  guantity  called 
the  thernal  conductance.  One  may  look  at  the  equation  as  an 
analogue  cf  the  basis  electrical  eguation  I  =  E/R  where  I  is 
eguivalent  tc  Q,  E  is  the  same  as  AT^  and  UA  is  1/fi.  This 
analogy  is  used  when  developing  the  expression  for  UA. 

In  order  to  account  for  the  effects  of  the  tube 
seeing  many  nearby  neighbors,  the  equation  is  modified  by 
tte   inclusion  of  a  factor,  F,  to  account  for  the  effects  of 
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the  surrcundir.g  tutes.   Equation  16  is  changed  to 


Q  =  u*A  F  AT 


m 


(17) 


where  F  is  defined  as  the  ratio  of  the  heat  flux  of  the  tube 
hurdle  tc  the  heat  flux  for  a  single  tube  in  a  stagnant  pool 
boiling  situation.  According  to  Ref.  13,  a  report  published 
by  workers  at  Heat  Transfer  Research,  Inc.  (HTBI)  of 
Alhambra,  California,  F  is  greater  than  cne.  In  figure  5, 
which  is  sketched  from  a  figure  in  their  report,  F  ranges 
frcu  U. 73  at  ATW  =  8°F  to  2.63  at  ATm  =  40°F.  Since  nc  data 
is  given  en  the  conditions  of  the  test,  F  is  assumed  to  be 
fixed  at  three. 

Ir  crder  to  find  AT^,  the  heat  exchanger 
ef f ectiveress ,  6  ,  is  determined.  The  effectiveness  is 
defined  as  the  ratio  of  the  actual  heat  transfer  rate  to  the 
maximum  ccssible  neat  transfer  rate  as  shown  by 


-mi 


(18) 


where 


C  =  (M  c>)nu(-j 


(19) 


CMjh  =  minimum  of  C^  or  Cc 


(20) 


In  the  bciler,  T^ jn  is  the  temperature  of  the  entering 
seawater,  ard  Tc in  is  the  saturation  temperature  cf  the 
wciking  fluid  at  the  pressure  at  the  mid  height  of  the 
bciler.  Beference  14,  by  Kays  and  Icndon,  contains 
eguations  for  finding  the  effectiveness  fcr  various  flow 
configurations.    In   addition  to  finding  Cmjn ,  Cmax  must  be 
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kncwn  in  crder  to  decide  which  equation  applies.  For  the 
case  of  ceiling,  C^cu  ^G  infinite;  therefore,  regardless  of 
th€  flew  pattern 


E  =  1  -  EXP(-iJIU) 


(21) 


wheie      NIC,      the      number      of      heat      transfer    units    cf    a    heat 
exchanger,    is    defined   as 


NTD    =    UA/C^jn 


{22) 


and 


Cmin    ~    Mryy     C 


1W    ^fc&w 


(23) 


AT^is    found    from 


£  = 


UA   ATry, 


rvur\ 


AT, 


(24) 


and  the  teirpeiature  of  the  sea  water  at  the   outlet   cf   the 
tute  is  deteinined  frcm 


i    = 


AIL 


(25) 


where 


ATo  ~  Ih,in"  Tc,in 


(26) 


Ihe   guantity   UA   is   determined   next.    Figure 
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illustrates  th€  nethod  used.  If  UA  is  pictured  as  the 
reciprocal  cf  resistance  then  the  flow  cf  energy  frcn  one 
fluid  tc  the  ether  can  be  pictured  as  a  series  electrical 
circuit  and  UA  is  determined  from 


UA  =  1/(R1  *  R2  ♦  H3  +  R4)  (27) 

Reference  15,  by  Holman,  shows  how  each  resistance  can  be 
found.  B1  is  the  thermal  resistance  due  tc  convecticn  on 
the  inside  cf  the  tubes,  the  sea  water  side. 


1 

R1  =  (28) 

7ThSw  dj  L 


R2  is  the  thermal  resistance  of  the  tube  wall. 


ln(d0/d-() 

R2  =  1 ! —  (29) 

2  7T  kw  L 


R3   is   the   thermal  resistance  attributed  tc  boiling  en  the 
outside  cf  the  tufces. 


R3  =  (30) 


R4  is  the  tberaal  resistance  due  to  any   other   resistances, 
such  as  fculing  and  corrosion  deposits. 

Ihe  heat  transfer  coefficient  on  the  seawater  side, 
h&w ,  is  determined  from  the  dimensionless  group  called  the 
Nusselt  number  defined  as  follows  for  flow  inside  of  tubes. 


^0 


Nu 


_  \w    ^i 


'SW 


(31) 


In  CEdcr  tc  iind  the  Nusselt  number,  three  distinct  regions 
cf  flow  are  allowed  for  in  the  problem:  laninar, 
transiticr,  and  turbulent  flow.  The  division  points  are  in 
terms  of  values  of  the  Reynolds  number,  Re  =  udp/p , 

laninar:         0  <  Re  <  2,CC0 
transition:  2,000  <  Re  <  10, COO 
turbulent:  10,000  <  Re  <  infinity 

In  the  laninar  region  the  Sieder-Tate  correlation  is  used. 


NuL  =  1.86  (fle  Pr  dj  /L  ) 


V3 


(32) 


In  the  turbulent  region  the   Eittus-Boelter   correlation   is 
used . 


KuTu  =  0.02  3  Re 


0.8   0.3 
Pr 


(33) 


In  the  transition  region  the  flew  is  unstable,  but  fron  data 
in  Ref.  14  it  appears  that  a  fairly  smooth  curve  drawr  from 
the  laminar  curve  to  the  turbulent  curve  adaguately 
represents  the  value  cf  the  Nusselt  number  in  the  transition 
regicn.  Ibis  is  done  by  assuming  the  form 


KuTr  =  y  NuL  +  (1  -  y)  NuTu 


(34) 


where 


y=aRe+tRe+cRe+d 


(35) 
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Th€  coefficients  in  the  polynomial  are  found  ty  requiring 

NuJr  =  NuL  at  Re  =  2,000 
NuTr  =  NuTu  at  Re  =  10,000 

and   that  the  first  derivatives  with  respect  to  the  Reynolds 
nunrer  match  at  Re  equal  to  2,000  and  10,000. 

The  I.I.  Mostinski  correlation,  from  Ref.  16,  for 
single  hcrizcrtal  tube  pool  boiling  is  used  to  find  the  heat 
transfer  coefficient  on  the  outside  of  the  tube.  HTRI's 
extensive  experimentation  led  Palen,  Yarden,  and  Tafccrek,  in 
fief.  13,  tc  conclude  that  Mostinski's  eguaticn  produced  the 
••mcst  consistent  results".   His  correlation  is 


C.69  0.7       0.17      1.2       1C 
hwj  =  C.CC658  pc    q    (1.8  r     +  4  r    +  10  r   )  (36) 

where 

r  =  F/Pc 

pc  =  critical  pressure  of  the  working  fluid  (psia) 

g   =  heat  flux  (BTU/hr-f t2) 

A  final  heat  flux  relation  needed  is  the  relation  fcr  the 
maximum  heat  flux  allowed.  Mostinski,s  ccrrelation  fcr  a 
single  tune  is  used. 


0.35       0.S 
q_   =  803  p„  r     (1  -  r)  (37) 

where  the  units  are  the  same  as  in  equation  (36).  From 
HTFI's  data,  under  unknown  conditions,  qKV,ax  for  the  rundle 
appears  tc  he  cne  third  of  the  qmQK  for  a  single  tube. 
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Ire  fewer  required  to  pump  the  salt  water  through 
the  tubes  defends  en  the  pressure  drop  across  the  tufce  bank. 
There  are  three  causes  of  the  pressure  drop.  The  flow 
experiences  a  loss  as  the  water  enters  into  and  exits  from 
the  tubes  in  addition  to  the  frictional  loss  along  the  tube. 
Ficm  Ref .  14,  the  pressure  drop  for  a  liquid  is  determined 
f  rem 

2    U  f  L 

AE  =  0.5  piw    uswv  (  __ ♦  Kc  -  Ke)  (38) 

di 

In  crder  tc  determine  f,  Ke  and  Kc,  the  flew  regime  must  be 
specified,  lie  regimes  are  a  little  different  than  those  for 
the  Nusselt  rumber. 

laminar:  0  <  Re  <  2,0C0 
transition:  2,000  <  Re  <  5,000 
turbulent:   5,000  <  Re  <  infinity 

The  Fanning  friction  factor,  f,  in  the  laminar  regime  is 

f  =  16/Re  (39) 

In  the  turtulent  regime,  the  Fanning  friction  factcr  is 
found  frcor  the  Blasius  equation  in  the  range  5,O0C  <  Be  < 
1CC,000 


0.25 
f  =  0.079/Re  (40) 

which   is   a   close   approximation  to  the  Karman  -  Nikuradse 
equation , 
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-0.5  0.5 

<4  f)      =  -0.6  +  0.87  ln(Re(4  f)    ) 


<41) 


which  is  valid  throughout  the  turbulent  region  (Re  >  5,000), 
assuming  smccth  called  tubes.  Either  eguaticn  (40)  cr  (41) 
may  be  used  with  similar  results  but  the  Elasius  is  letter 
suited  tc  this  analysis  because  the  Karman-Nikuradse 
eguation  mist  be  solved  by  iteration.  Eguaticn  (41) 
consumes  nore  computer  time  to  solve  and  can  cause  accuracy 
problems  if  not  enough  significant  digits  are  evaluated. 
The  fricticr.  factor  in  the  transition  region  is  determined 
by  fitting  c  least  sguares  curve  through  a  hand  drawn  curve 
from  the  end  cf  the  laminar  curve  to  the  beginning  of  the 
turbulent  curve. 

The  entrance  effect,  Kc,  and  the  exit  effect,  Ke, 
depend  on  the  Reynolds  number,  the  free  flow  to  frontal  area 
ratio,  <t>  ,  and,  in  the  laminar  region,  on  the  tube  length. 
For  the  boiler,  <t>  is  determined  by 


♦  - 


d;  N 


Ds   -  4  SA 


|42) 


Figure  S,  page  7-13,  in  Ref.  17  is  used  tc  find  Ke  aDd  Kc. 
Tfce  effect  cf  Reynolds  number  and  L/d  ratio  was  eliminated 
from  the  relationships  in  the  laminar  and  the  turbulent 
regimes  cf  the  flow,  but  the  Reynolds  number  was  included  in 
the  transition  region  because  of  the  large  variaticn  of  Kc 
with  Reynolds  number.  The  following  eguaticns  are  used  to 
find  the  entrance  effect  in  the  various  flow  regimes. 

laminar:  (using  the  curve  for  4  (L/d) Re  =  0.1) 


KcL  =  1.  -  0.4  <t) 


(43) 
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Surbulent:  (using  the  curve  for  Re  =  5,0CC) 

KcTu  =  0.52  -  0.4  <f>  (44) 

Iransiticr : 

KcTp  =  A  KcL  +  (1  -  A)  KcTu  (45) 

where 

A  =  e  He  +  f  (46) 

The  coefficients  e  and  f  are  determined  ty  requirirg  the 
transition  Kc  to  egual  the  laminar  Kc  at  Re  equal  to  2,000 
and  the  turtulent  Kc  at  fie  egual  to  5,000.  Ihe  exit  effect 
is  the  sane  icr  all  flow  regimes,  and  is  found  frcn  the 
following  eouation  that  approximates  the  4  (L/d) Re  =0.05 
curve  in  Bef .  17. 


2  3 

Ke  =  1.CCS3  -  2.5178<f>  +  1.1613(f)   -  0.  17677  <j>      (47) 


2 •  IhS    Cc nden ser 

lie  design  of  the  condenser  is  almost  the  sane  as 
the  design  cf  the  boiler.  It  is  a  fixed  tube  sheet, 
tuhe-in-stell  heat  exchanger.  The  condenser  does  net  reed  a 
separation  vclume  so  the  entire  shell  is  filled  with  tubes. 
In  crder  to  allow  the  entering  vapor  to  reach  all  cf  the 
tubes,   channels   through   the   tube   bank   are  needed.   The 
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allowance  for  this  is  accomplished  by  decreasing  the  shell 
dianeter  by  a  fixed  factor  when  calculating  the  number  of 
tubes.  The  calculation  of  the-  working  fluid  side  heat 
transfer  coefficient,  h  ,  is  different.  The  equation 
accounts  for  the  condensate  from  one  tube  falling  en  those 
belcw  it  ard  for  the  subcooling  of  the  condensate.  M.M. 
Chen's  equation,  as  developed  in  Eef.  18  and  shown  in  a  more 
sinple  fciir  in  Ref .  17,  is  used. 


^■"■[-■r)H[",.,t^i 


0.25 


(48) 


where 


h    =    BlU/ft2-hr-°F 


AT    =    Tsat     -    T^ 


(4S) 


H»       =    H,        +   —  CbAT 

tg         'g        8    p 


(50) 


n  =  number  of  tubes  in  a  column 


H'  is  the  enthalpy  corrected  for  subcooling  of  the 
ccrdensate.  Eguaticn  (48)  replaces  equation  (36)  in  the 
analysis  cf  the  condenser.  In  the  circular  shaped  tube 
tar.k,  n  is  fcurd  by  determining  the  number  cf  tubes  in  a 
column  of  average  height. 


Ds 
n  =  9.425  —  -  1 
Pi 


(51) 


The  definiticn  of  AT0  becomes 


k6 


±?o     =  l««t   -  I 


swjk 


(52) 


in  ccctrast  tc  equation  (26)   above 
temperatuie  is  found  from 


lhe   seawater   cutlet 


£  =   sw/out 


Lsw,»n 


AT, 


(53) 


which  replaces  equation  (25)  . 

In  cider  to  find  the  heat  transfer  ccefficient,  the 

temperature   cf   the   outside   tube   surface   is   needed  in 

equation  (49)  ;  therefore,  the  thermal  netwcrk  is  solved  for 
the  outside  kail  temperature  and  for  the  heat  flow  rate. 


Q  =  DA-  <TW  -  Tbulk  ) 


(54) 


where 


UA«  =  1/(R1  +  R2  +  R4) 


(55) 


Tb  ii  =  average  of  the  inlet  and 
the  outlet  seawater  temp. 

There  is  nc  factor,  F,  in  the  heat  transfer  equaticn  for 
ccrdensation  since  the  effect  of  neighboring  tubes  is 
acccunted  fcr  in  equation  (47)  . 


3.         3i.€_jrUJL£S 

Ihere    are    three    pumps    in   the   model.  Eoth 

circulaticn  pumps  for   the   boiler   and   the   condenser  are 

assumed   to  he  propeller  pumps  (axial  flow  pumps).   The  feed 
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pumps  are  assumed  tc  be  centrifugal  pumps.  The  differences 
in  the  characteristics  of  the  twc  types  of  pump  show  up  only 
in  the  calculaticn  of  capital  costs  of  the  pumps  ir.  the 
model. 

lie  chcice  of  the  type  of  pump  to  apply  is  made  on 
the  basis  cf  the  flow  rate,  GEM,  and  on  the  head  the  pump  is 
working  against.  Qualitatively,  the  propeller  pump  is  used 
where  the  flew  rates  are  high  and  the  head  is  lew.  The 
centrifugal  pumps  apply  where  the  flow  rate  is  comparatively 
lc*  and  the  head  pressure  is  high.  The  qualitative  picture 
is  placed  in  tetter  perspective  through  the  calculaticn  of 
th€  specific  speed  frcm 


0.5 
BEM   GEM 
Nsg  =  (56) 

0.75 
Head 


The  centrifugal  pump  is  used  when  500  <  Nsg  <  7,500  and  the 
propeller  puop  applies  when  7,500  <  Nsg  <  15,000,  according 
to  Bef.  1S.  The  choice  of  the  type  of  pump  is  checked  after 
the  optimi2aticn  process  has  been  completed  by  solving 
eguation  (55) .  If  the  value  of  the  specific  speed  is 
outside  cf  the  range,  500  <  Nsg  <  15,000,  the  maximum 
efficiency  cf  the  pump  drops  off  rapidly;  therefore,  few 
pumps  are  designed  tc  operate  outside  cf  this  range. 

Since  the  pumps  are  "black  boxes",  only  the  power 
consumed  by  the  pumps  is  calculated.  The  pewer  reguired  to 
drive  the  punp  is  found  from 


M  Ap 
P  =  __  ,57) 

fpurrtp 
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E.   COST  PFAMEKOFK 


As  discussed  at  the  start  of  this  report,  the  capital 
costs  of  the  major  components  are  the  only  costs  to  be 
considered  ir  this  mcdel.  The  boiler  and  the  condenser  are 
both  tu te-in-shell  heat  exchangers  assumed  to  te  cf  the 
fixed  tube  sheet  type.  The  only  difference  being  that  the 
shell  of  the  boiler  is  partly  filled  with  tubes.  Eoth 
seawater  circulation  pumps  are  assumed  to  be  propeller  type 
punps  and  nade  cf  the  same  material.  The  size  and  number 
per  heat  exchanger  shell  may  be  different  fcr  the  bciler  and 
the  ccndensei.  Ihe  number  of  working  fluid  feed  pimps, 
which  are  the  centrifugal  type  pumps,  are  assumed  tc  equal 
the  number  cf  bciler  shells. 

The  ccst  estimating  relationships  fcr  the  heat 
exchangers  ard  the  centrifugal  jumps  are  taken  from  Bef.  20, 
by  K.M.  Guthrie.  The  equation  fcr  the  cost  of  the  heat 
exchangers  is 


where  the  ccr.stants  are  calculated  from  a   leg-log   plct   of 
cost  versus  tctal  tube  surface  area  and 

Fj  =  0.6  for  a  fixed  tube  sheet 

Fw  =  C.C  for  pressure  correction  (p  <  150  psi) 

F^  =  naterial  factor  from  Table  1 

Ihe  =  ccst  index 

A^  =    d0  L  N 

fcr  a  relevart  range  cf  100  <  A   <  10,000  ft*. 

The   ccst   eguation   for   the  centrifugal  pumps  and  the 
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associated  meters  is 


0.721 
$f»  =  1'C13  Fm  Fo  I*>  CH  (59) 

wh€te  the  cocstants  are  found  from  a  straight  line 
approximation  cf  the  log-log  plot  cf  cost  versus  C  for  the 
section  cf  the  curve  from  C  =  30,000  to  C  =  300,000  for  an 
electric  motcr  driven  pump  and 

F0  =  1.0  fcr  suction  temperatures  <  250°F  (121°C) 

If.  =  cost  index 

F^  =  iraterial  factor  from  Table  2 

Cn  =  predict  of  the  flow  rate  in  GEM 

and  the  pressure  differential  in  psi 

for  a  relevant  range  of  30,000  <  CH  <  300,000.  The  cost 
cata  in  Guthrie's  article  is  given  for  a  time  base  of 
nid-1968. 

Reference  21,  also  by  Gutnrie,  gives  cost  data  for 
propeller  punps  and  centrifugal  pumps.  The  eguation  fcr  the 
prcpeller  runps  with  the  motors  is 


0.783 
$c^  =  5.767  Fm  ?p    Icp    Flow  (60) 

where   the  ccrstants  are  found  from  the  log-leg  plot  cf  cost 
versus  flew  rate  for  an  electric  motor  driven  pump  and 

Fp  =  1.0  fcr  a  suction  pressure  <  150  psi 

F^  =  naterial  factor  from  Table  3 

Ick  =  cost  index 

Flow  =  flew  rate  in  GPM 

for   a   relevant   range   of  1,000  <  Flow  <  100,000  GEK.   The 
tirce  base  fcr  the  prcpeller  pump  data   is   fcr   the   ecd   of 
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Table  1 

MATERIAL  FACTORS  FOR  THE  HEAT  EXCHANGERS 

Material 
Shell  Tube  Fm 

Carbcr  Steel      Carbon  Steel  1.00 

Carbcn  Steel      Brass  1.52 

Carbcn  Steel      Stainless  Steel       3.52 

Stainless  Steel   Stainless  Steel       4.50 
Carbcn  Steel      Monel  3.75 

flcnel  Monel  4.95 

Caibcn  Steel      Titanium  11.10 

litarium  Titanium  16.60 


Table  2 
KJTEEIAL  FACTORS  FOR  THE  CENTRIFUGAL  FEED  EUKP5 


Material 

Fm 

Cast  Iron 

1.00 

Bronze 

1.28 

Cast  Steel 

1.32 

Stainless  Steel 

1.93 

Monel 

3.23 

Hastelloy  C 

2.89 

Titanium 

8.98 
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Table  3 
MATEBIAL  FACTOBS  FOB  PBOPELLEE  PUMPS 

pumps 

Material  Em 

Cast  Iron  1 .00 

Cast  Steel  1.28 

Stainless  Steel  1 .64 
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1970.  Ike  costs  arrived  at  from  equations  (58),  (59),  and 
(6C)  are  in  dollars  per  component,  for  example,  dollars  per 
heat  exchanger  shell.  The  purpose  of  the  cost  incex  is 
alien  all  the  equipment  to  be  costed  on  a  common  basis.  The 
tine  difference  is  assumed  to  be  small  enough  so  that  the 
cost  index  is  set  to  1.00  for  all  the  components. 


E.   ASSEFEU  CP  THE  MATHEMATICAL  MCD2L 


The  engineering  and  the  ccst  relationships  are  linked 
together  in  the  model.  The  cost  eguation  is  the  function  to 
be  minimized  and  the  implicit  and  explicit  constraints  are 
the  engineering  eguations  that  restrict  the  feasible  region 
of  the  problem.  The  explicit  constraints  restrict  the 
variables,  tte  vector  X,  to  positive  values.  The  rest  of 
the  model  is  fcmulated  in  the  following  form 


mrimize 


f  (X) 


subject  to:   g  (X)  >  0,  1=1,2, ....  ,m 
i 

h  (X)  =  0,  i=m+1 , . . . . ,m+n 
i 


where  the  definitions  of  the  elements  of  X  are  shofcn  in 
Table  4.  Several  of  the  variables  must  have  only  integer 
values  if  the  system  were  actually  built,  but  for  the 
objectives  cf  this  aralysis  they  are  left  as  continuous  real 
variables.  For  example,  the  number  of  shells  may  ha\e  an 
integer  value  in  the  actual  system  but  this  requirement  can 
be  accomodated  by  solving  the  system  for  a  real  value  then 
integeri2ing  the  number  of  shells  and  solve  the  prcfclem  with 
the  next  higher  anc  the  next  lower  integer  number  of  shells. 
The  lower  cf  the  two  solutions  would  be  chosen.  The  final 
rescluticn   cf   the   correct    number    of   shells    would 
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Table  4 
DEJIKIIICN  CF  THE  ELEMENTS  OF  THE  VECTOR  X 


X  | 

[1) 

X  1 

[2) 

X 

[3) 

X 

[4) 

X  I 

[5) 

X 

[6) 

*  I 

7) 

X  1 

[6) 

X  I 

(9) 

X  I 

[K) 

X 

[11) 

X 

[12) 

X 

[13) 

velocity  of  seawater  in  the  boiler 
«•     »     »     "    "   condenser 

cutlet  pressure  in  the  boiler 

inlet     "     "    "   condenser 

tube  length  in  the  boiler 
11     "     ••    "   condenser 
mass  flow  rate  of  the  working  fluid 

diameter  of  the  boiler  shell 
"     "    "   condenser  shell 

outside  dia.  of  a  tube  in  the  boiler 
ii     ii    ii   ii   it    ii    ii   condenser 

number  of  bciler  shells 
"    "   condenser  shells 
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prcbably  defend  on  other  than  thermodynamic  and  cost 
considerations,  such  as  a  requirement  that  ere  extra  shell 
be  included  to  allow  for  one  shell  being  down  for 
maintenance.  Obviously,  the  number  of  tubes  in  a  shell  must 
be  an  integer  value  but  one  extra  tube  cut  of  several 
thousand  is  a  nincr  matter.  In  addition,  there  may  be  a 
requirement  fcr  extra  tubes  to  allow  for  a  percentage  tc  be 
plugged  when  they  develop  leaks  so  that  the  entire  shell  is 
net  put  cut  cf  ccmirission  by  a  minor  leak  in  one  tube. 

1«   Objective  Function  and  the  Constraints 

The  objective  function  is  the  sum  of  the  costs  of 
the  boiler  arc  the  condenser,  and  the  pumps.  The  model 
assumes  ttat  there  is  one  circulating  pump  and  one  feed  pump 
for  each  shell  in  the  bciler  and  that  there  is  one 
circulating  pump  for  each  shell  in  the  condenser.  Since 
there  is  nc  cost  analysis  for  the  turbine,  nc  assumpticn  is 
needed  abcut  the  number  of  turbines  in  the  system. 

There  are  nineteen  constraints  in  the  present 
prcgram,  cne  cf  which  is  optional.  The  first  three  cencern 
the  saturaticn  pressures  in  the  boiler  and  the  condenser. 
The  next  twelve  place  restrictions  on  various  dimensiens  of 
the  bciler  anc  the  condenser.  The  last  three  require  that 
scne  intermediate  calculations  stay  positive. 

Using  the  objective  function  as  a  constraint  has 
proved  advantageous  in  cases  where  the  problen  was 
nen-cenvex  or  where  the  program  had  difficulty  in  finding  a 
solution.  It  is  used  by  setting  a  boundary  en  the  objective 
function  so  that  the  constraint  is  either  infeasitle  or 
slightly  feasible  at  the  starting  point.  The  proper  choice 
must  be  fcund  by  trial  and  error.  One  caution  to  note  is 
when   equality   constraints   are   involved.   If  the  starting 
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point  is  neat  the  optimum,  the  program  tries  to  move  away 
frcm  the  boundary  and  the  equality  constraints  may  not  be 
satisfied  wten  the  program  finishes. 

The  need  for  the  ether  eighteen  ineguality 
constraints  is  net  obvious  when  the  problem  is  first 
formulated.  The  thermodynamics  of  the  problem  reguire,  for 
example,  that  the  propane  saturation  temperature  be  less 
than  seavater  temperature  in  crder  for  the  toiler  to  absorb 
energy,  tut  the  program  has  no  way  of  knowing  that,  unless 
it  it  specified.  If  the  program  can  improve  the  cost  and 
satisfy  the  ether  constraints  by  reversing  the  relationship 
it  will  dc  so.  This  inversion  may  take  place  only  for  a 
short  tine  during  the  solution  of  the  problem,  but,  if  the 
eguations  fail  under  this  condition,  the  constraint  en  the 
temperatures  in  the  bciler  is  necessary.  Most  of  the  other 
constraints  are  needed  for  the  same  reason. 

Seme  constraints  are  not  added  to  the  problem  since 
they  are  net  binding  in  the  solution.  There  are  two  such 
constraints.  The  first  is  the  segmental  area  reguired  to 
allow  separation  of  the  liguid  droplets  from  the  vapcr.  The 
segmental  area  set  aside  by  the  H0%  rule-cf-thumb  is  more 
thar  that  reguired  by  the  consideration  of  the  physical 
properties  of  prcpane.  The  ether  constraint  is  the 
restriction  en  the  maximum  heat  flux  by  Mostinski's 
correlation,  eguation  (37).  The  constraints  dc  not  tecome 
binding  in  tfce  solution. 

The  three  eguality  constraints  contain  the 
calculations  of  the  properties  of  the  various  components 
modeled  in  this  analysis.  The  first  cne  concerns  the 
bciler.  It  states  that  the  actual  energy  absorption  rate  of 
the  bciler  shells  must  be  egual  tc  the  energy  absorption 
rate  reguired  by  the  thermodynamic  cycle.  Similarly,  the 
second  eguality  constraint  states  that  the  heat  rejected   by 
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the  condenser  must  equal  that  required  by  the  thermod vnamic 
cycle  for  the  condensation  of  the  vapor  from  state  point  5 
back  to  state  pcint  1  as  shown  in  Fig.  3.  The  third 
ccnstraint  calculates  the  pcwer  required  by  the  various 
punps  and  the  pcwer  extracted  from  the  turbine  to  run  the 
system's  punps  and  to  produce  electrical  power  for 
ccnsumpti.cn.  This  ccnstraint  requires  that  the  power  cut  cf 
the  turbine  equal  the  power  used  by  the  jumps  and  the 
generator.  The  net  useful  power  cut  of  the  systen  is  a 
design  cfc^ective  and  is  specified  as  25MW  herein. 

2«   Ibe  Optimization  Method 

Since  the  arrival  of  the  electronic  computer,  a 
great  nany  nonlinear  programming  techniques  have  teen 
developed  frcm  early  theories  and  more  current  research. 
The  ccffputei  program,  called  the  Sequential  Unconstrained 
Minimization  Technique  (SUMT) ,  developed  by  Mylander, 
Helens,  and  McCormick  in  Bef.  22  is  used  to  optimi2e  the 
ffcdel.  Their  prcqram  implemented  much  cf  the  theory 
contained  ir  Ref.  23,  by  Fiacco  and  McCormick,  about 
ncrlinear  programming  using  unconstrained  minini2ation 
techniques.  SUMT-Version  4  is  chesen  because  of  its  ability 
to  solve  a  wide  variety  of  problems.  The  theoretical 
requirements  for  SUMT  to  be  guaranteed  to  find  a  local 
minimum  is  that  there  exists  seme  point  which  satisfies  the 
ineguality  constraints  and  that  there  not  be  a  local  minimum 
at  points  where  X  goes  to  infinity.  If  certain  convexity 
conditions  are  satisfied,  then  the  local  ninimum  found  is  a 
glcbal  minimun. 

SIMT  uses  a  penalty  function  to  transform  the 
censtraired  cptimizaticn  problem 
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nicimi2e     f  (X) 

subject  tc:   g  (X)  >  0,  i=1r2,  . . . .  ,m 
i 

h  (X)  =  0,  i=m+1 , . . . . ,m+n 
i 


intc  a  sequence  of  unconstrained  problems  of  the  form 


m  m  +  n 

E(X,r)  =f(I)  -r.Ilng  (X)  ♦  1/r,  I   <h.  (X))  2    (61) 

i=1     1         i=m+1   1 

where  E(X,r)  becomes  the  function  to  be  itinimized.  The 
inequality  constraints,  g  (X)  ,  cause  a  large  penalty  to  occur 
as  a  boundary  is  approached.  The  equality  constraints  add  a 
large  penalty  whenever  h  (X)  departs  from  zero  in  either 
direction.  The  size  of  the  penalties  due  tc  g  (X)  and  h (X) 
is  controlled  by  the  arbitrary  parameter  "r".  The  algorithm 
solves  the  problem  by  solving  a  series  of  sutproblems.  For 
r1 >r^>r3  . .  .>0,  SUMT  minimizes  P(X(rh),rn)  in  each  sukproblem 
in  crder  to  find  the  minimum  of  f  (X)  as  "r"  tends  to  zero. 
When  a  minimum  of  a  subproblem  is  found,  "r"  is  reduced  by  a 
factor  assigned  by  the  user  and  the  new  subprcblem  is 
solved.  The  process  continues  until  the  convergence 
criteria  chosen  by  the  user  is  satisfied.  Figure  7 
illustrates  schematically  the  results  of  the  seguence  of 
sutprcblens  and  the  qualitative  shape  of  the  penalty 
function.  Three  methods  of  minimization  are  included  in  the 
program.  The  user  may  choose  either  the  generalized 
Newton-Eaphscn,  steepest  descent,  or  McCormick  modification 
of  the  Fletcher-Powell  method. 

Ihere  are  several  other  options  that  control  the 
operation  of  SUMT.  The  user  may  cheese  the  initial 
penalty,  1,  cr,  if  the  problem  satisfies  certain  conditions, 
the  user  nay  elect  to  have  SUMT  compute  the  initial  "r". 
There  are  three  choices  for  the  convergence  criterion  of  the 
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FIGURE  7.  PROBLEM  SOLUTION  SEQUENCE 
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problem.  lh€  criterion  used  depends  on  the  inequality 
penalty  tern  teing  close  enough  to  zero.  Another  option  is 
the  choice  of  the  calculation  of  the  subproblem  stepping 
criterion  .  Two  of  the  choices  compute  the  inverse  of  the 
second  partial  matrix  and  the  third  choice  computes  only  the 
gradient  of  the  penalty  function  (  VxP(X,r)).  One  of  the 
most  useful  options  is  extrapolation.  After  two  or  three 
suhproblems  have  been  solved  SUMT  is  able  to  make  a  first  or 
second  crcer  extrapolation  to  the  starting  point  for  the 
next  sutprcblem.  This  can  speed  up  the  solution  of  the 
problem  measureably. 

Ike  characteristics  of  the  problem  limit  the  choices 
for  each  option.  Ihe  penalty  factor  couldn't  be  set  bj  SUMT 
due  to  the  presence  of  equality  constraints  and  because  the 
iritial  3!  is  usually  not  close  to  any  boundary.  The  stepping 
criterion  for  the  problem  is  that  the  inequality  penalty 
tern  beccme  less  than  some  value,  usually  10~7. 
Extrapolation  is  used,  but  it  makes  little  difference  in  the 
answer  or  to  the  solution  time  whether  a  twe  point  or  three 
point  estinate  is  made.  The  subproblem  convergence 
criterion  is  the  magnitude  of  the  gradient  of  I  (X,r) 
becoming  less   than  some  small  number. 

Ike  choice  cf  starting  penalty,  r,  affects 
significantly  the  solution  time  and  even  whether  a  solution 
is  found  it  seme  problems  according  to  Eef.  7.  If  the 
initial  "r"  is  toe  small,  the  subproblem  optimum  may  he  the 
prcrlem's  optimum  which  may  be  hard  to  find  or  require  more 
calculations  than  if  a  much  larger  "r"  is  chosen.  If  "r"  is 
toe  large,  the  first  few  optima  lie  near  the  center  of  the 
feasible  region  and  are  unrelated  to  the  optimum  cf  the 
objective  function. 

Scaling  of  the  variables,  constraints,  and  the 
objective  function  is   necessary.    The   variables   must   be 
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within  tkc  cr  three  orders  of  magnitude  of  each  other,  for 
example,  the  working  fluid  flow  rate  is  of  the  order  of  107 
whereas  the  tube  diameters  are  of  the  order  10°.  If  the 
finite  difference  interval  is  10— ♦,  the  partial  of  the 
equations  with  respect  to  the  mass  flow  rate  may  be  lost  due 
tc  the  rcundcff  error  providing  the  only  significant 
figures.  with  the  constraints  and  the  objective  function, 
unsealing  may  result  in  one  term  in  the  penalty  furction 
being  deninant,  and  the  algorithm  finds  that  most  cf  the 
decrease  in  the  penalty  function  takes  place  in  the  dcninate 
term.  The  result  is  that  SUMT  wcrks  hard  tc  decrease  that 
term  and  does  not  see  the  other  terms  resulting  in  a 
meaningless  scluticn. 

Ir  making  first  and  second  order  estimates  cf  the 
scluticn,  SLMT  may  generate  negative  numbers  in  the  X 
vector.  When  it  evaluates  the  constraints  and  objective 
furcticn  for  this  point,  the  equations  fail  and  may  cause 
the  program  tc  fail  depending  en  the  severity  of  the  failure 
and  how  many  times  the  error  has  cccured.  This  problem  is 
circumvented  by  testing  the  vector  upon  entry  intc  the 
subroutine  where  the  objective  function  and  the  constraints 
are  evaluated  and  setting  the  constraint  being  evaluated 
equal  to  seme  negative  constant. 

Ihe  main  program  of  SUMT  is  altered  so  that  the 
prcblem  can  be  restarted  at  the  current  solution  without 
reading  in  a  new  problem.  This  is  necessary  because  the 
prcgram  would  give  all  the  indications  cf  having  fcund  a 
solution,  but  if  the  problem  is  restared  it  cculd  ofter  find 
a  lew  value  for  the  objective  function.  The  new  soluticn  is 
compared  to  the  previous  solution  and  if  the  two  solutions 
are  net  clcse  enough  the  problem  is  restarted  after  naking 
seme  changes  to  particular  parameters,  optiens,  or 
constraints  depending  on  the  objective  the  user  has  in  mind. 
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Ill 


RESULTS  AND  CONCLUSIONS 


One  cf  the  majcr  tasks  in  the  development  of  a  nonlinear 
model  is  the  linking  of  the  model  with  the  optimization 
program  in  order  to  analyze  the  model.  The  selection  cf  the 
prefer  ncclinear  programming  technique  is  essential  to  the 
solution.  Several  techniques  require  that  the  objective 
function  and  the  constraints  be  in  a  particular  fcrm.  For 
example,  Geometric  Programming  developed  by  Euffin,  Peterson 
and  Zenei,  in  Ref.  24,  reguires  that  the  eguatiens  be 
pcsyncmials  *hich  is  not  true  of  the  model  developed  ir  this 
paper.  SUM!  was  chosen  because  cf  its  wide  applicaticns  to 
problems  cf  any  tyre.  It  only  reguires  that  there  exist 
some  feasible  point  and  that  local  minima  not  exist  at 
infinity.  It  has  successfully  solved  many  nenconvex 
problems . 

Twc  versiens  of  the  analytical  model  are  presented.  The 
first  is  called  the  iterative  version  because  the 
calculaticn  cf  the  heat  absorption  rate  and  heat  rejection 
rate  reguire  iterating  on  some  quantity  to  find  them. 
Mostinski's  correlation,  eguation  (36),  contains  the  heat 
flux  which  is  an  unknown.  It  is  not  possible  to  solve  the 
system  cf  equations  explicitly  for  the  heat  flux;  therefore, 
a  successive  approximation  technique  is  used.  The 
calculaticn  cf  the  heat  rejected  in  the  condenser  involves 
the  same  piocess  with  the  temperature  of  the  outside  tube 
wall  as  the  iterating  variable.  SUMT  is  urable  to  arrive  at 
a  consistent  solution  in  the  time  allowed,  up  to  90  minutes 
of  computer  time.  The  solutions  are  quite  different  for 
each  try  vith  nc  pattern  evident  about  the  directicn  the 
scluticn  night  lie. 
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There  are  several  possible  explanations  fci  this 
behavior.  Eerhaps  the  behavior  of  this  version  cf  the 
problem  is  such  that  SUMT  is  not  the  best  choice  fcr  the 
optimizing  algorithm  or  possibly  the  prchlera  cannct  be 
solved  at  all  ir  this  form.  Another  source  of  problems 
could  be  the  numerical  differentiation  procedure.  Numerical 
differentiation  is  subject  to  severe  truncation  and 
rcund-cfi  error.  The  subproblem  always  terminates  when  SUMT 
is  unable  to  reduce  the  penalty  function  during  the  next 
iteration,  frcm  one  estimate  of  X  to  the  next  estimate. 
Eemember  the  convergence  criterion,  never  satisfied  for  this 
version  cf  the  problem,  is  the  magnitude  cf  VxE(X,r) 
becoming  less  than  some  small  number.  No  conclusive  results 
could  re  obtained  with  this  version  so  the  second  version  is 
used  fcr  tbe  present  analysis. 

In  this  second  version  the  heat  transfer  coefficierts  on 
tbe  Ercpane  side  of  the  tube  are  fixed.  This  enables  the 
eguality  constraints  for  the  boiler  and  the  condenser  to  be 
calculated  to  the  full  accuracy  of  the  computer  without 
increasing  the  time  reguired.  The  correlations,  eguation 
(36)  and  (48) ,  are  used  for  comparison  with  the  fixed  value. 
If  the  calculated  value  and  the  fixed  value  are  within  20  to 
25%  of  one  another,  the  solution  would  be  considered  to  be 
gccd  since  this  is  the  tolerance  claimed  for  the 
correlations . 

SuMT  seemed  to  be  able  to  solve  the  new  problem  because 
a  series  cf  seven  solutions  have  been  produced  wheie  the 
values  of  tbe  elements  of  X  and  cf  the  objective  function 
vary  no  more  than  3.7%  and  1.2X,  respectively,  about  the 
mean  when  different  starting  points  and  different  SUMT 
control  parameters  are  selected.  See  Table  6  for  a  summary 
of  the  solutions  obtained  from  the  assumed  data,  as  shewn  in 
Table  5.  The  design  implications  of  these  results  are  guite 
reasonable.   The  range  shows  the  spread   in   the   solutions. 
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Table  5 
PROBLEM  INPUT  DATA 


Fouling  Thermal  Besistance 
Ecil€r  and  Condenser 

Elevation  cf  Boiler 
Aleve  Condenser 

Channel  Allowance 

Efficiencies  for  the  Turbine 

and  the  Pumps 
Seavater  Temperatures  -  Hot 

-  Cold 

Net  Output  Pcwer 

Eundle  Pactor  (F) 
Ccst  Irdexes  (I) 
lute  Hall  Thickness 

Mir.  Distance  Between  Tufces 

Heat  Transfer  Coefficients 
fcr  Eciler  and  Condenser 


0.0  hr-cF/BT(J 
(0.0  °C/kW) 
15  ft 
(4.57  ib) 
0.83  ft 
(0.25  m) 
85* 

75°F 
(23.9°C) 
40°F 
(4.4°C) 

8.532x107  BTU/hr 
(25  MW) 
3.0 

1 .00  fcr  all  cemp 
0.028  in 
(0.071  cm) 
0.5  in 
(1.27  cm) 

145  BlC/ft2-hr-°F 
(0.822  kK/m2-oC) 


Materials 


Eciler  and  Ccndenser 

Feed  Pumps 
Circulation  Pumps 


Carbon  Steel  Shells 
Brass  Tubes 
Bronze 
Stainless  Steel 
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EXAMPLE  PEOBLEM  SOLUTION  SUMMABY 


Variable 

8anc[e 

Mean 

Var 

rJn  its 

Bciler : 

Velocity 

3. 323-3.347 

3.33  5 

0.7 

ft/sec 

(1.013-1.020) 

(1 .016) 

( a/  s  e  c ) 

El€S£Uie 

ajprox  0 

118.2 
(815.0) 

o.c 

lb/in2 
(fcp) 

lute  length 

39.  16-40  .00 

39.71 

2.  1 

f + 

A.  — 

(11.91-12.19) 

(12.10) 

(at) 

Shell  dia. 

24.12-25.00 

24.9  1 

3.5 

ft 

(7. 352-7.620) 

(7.59  3) 

(m) 

lute  dia. 

0. 756-0.761 

0.758 

0.7 

in 

(1. 919-1.934) 

(1.926) 

(cm) 

#  of  shells 

7.  141-7.409 

7.228 

3.7 

- 

Ccr  denser : 

Velocity 

2.793-2.823 

2.808 

1.  1 

ft/sec 

(0.851-0.860) 

(0.856) 

(ra/sec) 

Eressure 

91  .53-91.60 

91.56 

0.1 

lb/in2 

(631.  1-631.6) 

(631.3) 

(kP) 

lute  length 

39.35-40  .00 

39.78 

1.6 

ft 

(11.99-12.19) 

(12.12) 

(m) 

Shell  dia. 

44.69-45.00 

44.89 

0.7 

ft 

(13.62-13.72) 

(13.68) 

(m) 

lute  dia. 

0.7632-0.7686 

0.766  0 

0.7 

in 

(1.938-1. 952) 

(  I. 946) 

(cm) 

4  of  shells 

3.297-3.401 

3.331 

3.  1 

- 

Mass  flcv  rate 

22.  12-22.21 

22.  17 

0.4 

10*ltm/hr 

(10.03-10.07) 

(10.06) 

(10*kg/hr) 

System  ccst 

7.576-7.673 

7.607 

1  .2 

106  $ 
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Th€  arithmetic  mean  of  the  seven  solutions  is  calculated. 
The  difference  between  the  maximum  and  the  minimum  values  is 
shewn  as  a  percent  of  the  mean.  The  percent  variaticn  gives 
sc&e  idea  cf  how  well  the  solution  was  determined.  Note 
that  the  tube  length  and  shell  diameter  constraints  are 
binding  for  roth  the  boiler  and  the  condenser.  The 
sutprcblem  still  terminates  only  when  the  penalty  furction 
does  not  decrease  from  one  iteration  to  the  next.  The 
magnitude  cf  ^F(X#r)  is  not  becoming  small  enough  to 
satisfy  tre  convergence  criterion.  The  value  never  beccmes 
less  thar  C.C2  and  is  usually  greater  than  0.5  vhen  a 
sutprcblen  is  stopped  by  the  program. 

Useful  results  are  obtainable,  however.  The  engineer 
can  learr  much  about  his  problem  even  from  a  single  point. 
Fcr  exanple,  frcm  the  cost  breakdown  in  Figure  8,  the 
portion  cf  tte  costs  attributable  to  each  component  can  be 
seen.  The  heat  exchangers,  of  course,  dominate  with  18%  of 
the  cost.  The  circulation  pumps  are  about  21$  of  the  cost 
with  the  feed  pumps  being  T.2X  of  the  total  cost.  Since 
the  efficiency  assumed  for  the  circulation  pumps  is  prcbably 
toe  high,  these  pumps  would  be  even  mere  ccstly  in  teims  of 
energy  and  meney.  This  could  indicate  the  priorities  that 
should  fce  set  on  the  design  efforts  for  improvements  cf  the 
vaiicus  cempenents  with  the  highest  priority  beinc  the 
highest  percentage  cost  component. 

The  syst€B»s  fewer  distribution  is  shown  in  Table  7. 
The  sigrificant  features  are  the  large  energy  absorption 
rate  for  the  snail  useful  output  power  and  the  parasitic 
pumping  power.  The  lower  the  efficiency  the  greater  must  be 
the  power  absorbed  for  the  same  output.  This  is  to  be 
expected  for  this  kind  of  system.  The  parasitic  pumping 
pewer  (1251  of  the  turbine  output)  is  about  twice  what  it  is 
for  a  fessile  fuel  plant. 
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FEED 
PUMPS 
$94,100 
1.2% 


BOILER 

CIRCULATION 

PUMPS 

$654,000 

8.6% 


CONDENSER 

CIRCULATION 

PUMPS 

$978,000 

12.91 


FIGURE  8.  COST  BREAKDOWN 
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Table  7 
SOLUTION  POWEE  EUCGET 


F€at  Absorption 
Turbine  Shaft  Power 
fe€d  Eumps 

Eli  Circulation  Pumps 
Ccfld  Circulation  Pumps 

E€at  f€J€Ct€d 

Ket  Iseful  Power 


ETU^hr 

MWatts 

3.56x1C« 

1042.8 

9.69x107 

28.4 

4.49x10* 

1.3 

2.58x10* 

0.8 

4.51x10* 

1  .3 

3.47x109 

1015.7 

8.53x107 

25.0 
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Cnce  the  mcdel  has  been  optimized  with  the  number  of 
shells  as  real  variables,  the  problem  can  be  re-solved  twice 
by  fixing  the  shells  at  the  next  higher  and  the  next  lower 
integer  ruirfcer.  For  the  higher  number  of  shells,  the  cost 
increased  slightly  as  expected,  from  7.607  million  dcllars 
average  to  7.983  million  dcllars.  The  problem  with  the 
lcwer  integer  number  of  shells  could  not  be  solved  due  to 
the  constraints  en  the  tube  length  and  the  diameter  cf  the 
shells  being  banding  before  the  shell  numbers  were 
integerized.  The  program  is  unable  to  find  a  new  solution 
that  satisfies  the  equality  constraints. 

Table  8  shows  some  miscellaneous  data  from  the  mcdel 
that  deserves  comment.  The  overall  heat  transfer 
coefficient  cf  the  bciler  is  approximately  three  tines  that 
for  the  cendenser  which  is  the  major  reason  for  the  lcwer 
ccst  of  the  boiler  relative  to  the  condenser.  However, 
remember  that  the  bundle  factor,  egual  to  3,  is  a  part  of 
the  heat  transfer  coefficient.  More  data  must  be  gotten 
abcut  this  factor  in  crder  to  improve  the  reliance  or.  any 
sclution.  There  is  a  significant  difference  between  the 
Prcpane  beat  transfer  coefficients  that  were  assumed  aid  the 
theoretical  values  calculated  from  the  correlations.  Other 
values  shculd  be  tried  until  the  fixed  value  is  withir  the 
20  5*  tolerance  on  the  theoretical  value.  The  flew  of 
seawater  tbrcugh  the  tubes  in  both  the  bciler  anc  the 
ccrdensei  is  turbulent.  Apparently,  the  resultinc  cost 
savings  due  to  the  improved  heat  transfer  rate  through  a 
smaller  area  is  more  than  the  increased  costs  of  the  higher 
pumping  pewer  needed.  It  seems  likely  that  the  solution 
obtained  is  globally  optimal;  several  different  starting 
pcints  were  tried.  The  problem  has  been  started  in  both  the 
laninar  and  the  turbulent  flow  regions  with  no  change  in  the 
results. 
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Table  8 
ADDITIONAL  DATA  FROM  SAMEIF 
PROBLEM  SOLUTION 

Boiler    Condenser     Units 


Working  fluid  Heat  Transfer  Coefficients 
Jssuraed 


145. 

145. 

BTU/hr-f t2-°F 

(0.822) 

(0.822) 

(kW/m2-cC) 

227. 

241. 

(1.29) 

(1.37) 

Seawater  Cutlet  Temp 
Shell  Saturation  Temp 


1 heoretical 


Surface  £rea  per  Sfcell   2.39x10s  1.27x106     ft2 

(2.22x10*)  (1.18x105)   (m2) 

72.3  43.8       °F 

(22.4)  (6.55)     (°C) 

66.7  49.9       °F 

(19.3)  (9.94)     (°C) 

Seawater  Pass  Plow  Rate   4.39x108  9.11x109  lbm/hr 

(1.99x108)  (4.13x109)  (kg/hr) 

6.67  7.85       °F 

(3.71)  (4.36)     (OC) 

310.6  105.0  ETU/hr-ft2-°F 

(1.76)  (0.595)  (kW/m2-°C) 

2070.0  824.9  ETU/ft2-hr 

(6.528)  (2.602)  (kW/m2) 

19,800.  9,960. 


Mean  Temp  Difference 

Overall  Feat  transfer 

Coef f ici  ert 
Heet  Flux 


Seawater  Eeyrclds  ro 


70 


IV 


RECOMMENDATIONS 


The  present  model  has  a  great  deal  mere  infomaticn  to 
offer,  but  scire  improvements  in  the  model  shculd  be  mace  to 
improve  the  ability  of  a  nonlinear  program  to  scl\ie  the 
prctlem.  The  first  recommendation  would  be  to  studj  the 
eguality  constraints  and  place  inequality  constraints  en  any 
vaiiable  cr  grouping  of  variables  if  they  rcth  can  tecome 
negative  ace  also  must  have  their  logarithms  computed. 
Examples,  already  included  in  the  program,  are  the  numter  of 
tufes  in  tie  shells.  The  heat  transfer  coefficients  fcr  the 
working  fluid  side  in  the  boiler  and  the  condenser  cculd  be 
included  as  equality  constraints.  This  can  be  dene  by 
adding  the  heat  transfer  coefficients  to  the  vector  X  and 
reguiring  that  the  variables  be  equal  to  the  coefficierts  as 
computed  ty  the  appropriate  correlation.  The  same  cculd  be 
dcre  with  Ealen  and  Small's  property-dependent  calculation 
of  the  separation  area  by  equation  (12)  and  (13).  Frcm  an 
overall  standpoint,  the  results  seem  to  indicate  either  the 
problem  is  tec  difficult  to  be  solved  efficiently  ir  its 
present  fern  or  that  another  nonlinear  program,  other  than 
SUMI,  might  re  better  suited  to  the  problem.  Numerical 
differentiation  of  the  constraints  in  this  problem  do  not 
seem  to  help  the  program  to  find  the  solution  since  the 
magnitude  of  ^LP(X,r)  is  net  near  zero  when  a  sutproblem 
steps. 

Once  the  optimization  scheme  is  functioning  properly, 
the  model  can  te  used  to  study  the  effects  of  fouling  and 
ccrrcsicr  and  the  costs  of  their  prevention.  Many  proposals 
have  been  made  that  must  be  analyzed  in  the  context  of  the 
total  sjsten.  Seme  such  as  Amertap,  where  plastic  foam 
balls  are   pumped   through   the   tubes   to   clean   the   heat 
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transfer  surfaces,  add  a  cost  penalty  free  the  capital  and 
operating  ccst  of  the  Amertap  system  and  frcm  the  increased 
pumping  pewer  reguired.  Coatings  add  the  capital  ccsts  of 
the  coatirg  and  the  added  capital  cost  of  increased  heat 
transfer  area  reguired  because  of  the  increased  thermal 
resistance  cf  the  coating. 

A  sensitivity  analysis  can  be  performed  en  the  model  to 
test  the  assumptions  and  the  various  constants.  The 
sensitivity  tc  changes  in  these  factors  can  indicate  how 
much  one  shculd  be  willing  to  pay  for  the  improvement 
sought.  High  sensitivity  indicates  areas  where  the  greatest 
effort  fcr  improvement  should  be  desired  or  the  priorities 
shculd  be  set  en  research  into  the  unknown  areas  cf  the 
problem.  The  trade-off  between  cost  and  performance  can  be 
studied  by  changing  the  output  power  and  looking  fcr  the 
lowest  ccst  per  kilowatt-hour,  which  is  the  measure  by  which 
ccrventicnal  pewer  plants  are  compared.  One  could  lcok  into 
the  cost  and  the  benefits  of  increasing  the  depth  frcm  vihich 
the  cold  water  is  pumped  versus  the  improved  temperature 
differential  possible  from  usirg  the  colder  water. 

For  the  capital  cost  of  the  system  to  be  complete,  the 
turbine  design  must  be  included  and  ccst  estimates  made. 
The  turbine  will  be  unique  due  to  the  very  low  pressure 
potential  drop  it  must  use.  Some  work  has  been  dore  by 
ether  grcups  working  on  OT£C,  Bef.  25  and  26.  Reference  27, 
by  Bcbert  1.  Bartlett,  can  provide  costing  and  design 
infermatien  although  he  deals  only  with  steam  turbines.  The 
remaining  part  of  the  system  that  is  desireable  for 
inclusicc  ir.  the  analysis  presented  herein  is  the  cold 
seawater  intake  pipe.  It  will  be  on  the  order  of  1,0C0  ft 
(3C4.6m)  lcr.g  and  100  ft  (30.4  m)  in  diameter. 
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